Distribution of sizes of erased loops for loop-erased random walks 
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We study the distribution of sizes of erased loops for loop-erased random walks on regular and fractal lattices. We show 
that for arbitrary graphs the probability P{1) of generating a loop of perimeter I is expressible in terms of the probability Pst{l) 
of forming a loop of perimeter / when a bond is added to a random spanning tree on the same graph by the simple relation 
P{1) — Pst(l)/l. On d-dimensional hypercubical lattices, P{V) varies as l^'^ for large /, where cr = 1 + 2/z for 1 < d < 4, where 
z is the fractal dimension of the loop-erased walks on the graph. On recursively constructed fractals with d < 2 this relation is 
modified to cr = 1 -I- 2d/{dz), where d is the hausdorff and d is the spectral dimension of the fractal. 

dimensions the latter is known using the equivalence of 
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The loop-erased random walk (LERW) is a simpler 
variant of the well-known self-avoiding walk (SAW) prob- 
lem, which retains the no self-intersection property of 
SAWs, but is closer to the ordinary random walk prob- 
lem. In this paper, we study the distribution of sizes of 
erased loops for LERWs on arbitrary graphs and relate it 
to the distribution of sizes of loops formed when a bond is 
added to a random spanning tree on the same graph. In 
the themodynamic limit this distribution has a power law 
tail, and we express the exponent in terms of the fractal 
dimension of the chemical paths on spanning trees. We 
show also how this relation is modified for deterministic 
fractals. 

In spite of the fact that the LERW model is some- 
what more tractable analytically than the SAW prob- 
lem, the number of papers devoted to this problem has 
remained rather small. The model was defined by Lawler 
, who called it the loop-erased self-avoiding walk (LE- 
SAW). This terminology is somewhat inappropriate and 
we prefer to use the term LERW in this paper 1^ . Lawler 
showed that for space dimensionality d > 4, the large 
length scale properties of the LERWs are same as those 
of simple random walks. Thus if r„ is the end-to-end 
distance of a n step LERW and we define the exponent ly 



such that < >^ 



then has the value 1/2 in d > 4. 



For d < 4, Lawler derived the rigorous bound that i/ > 
the flory value 3/(d+ 2) for SAW |]. From numerical 
simulations, Guttman and Bursill Q obtained the values 
V = 0.800±0.003 in two dimensions and v = 0.616±0.004 
in three dimensions. The corresponding values for SAWS 
are v ~ 0.750 and ly — 0.59 ± 0.004 in two and three di- 
mensions respectively. This shows that LERW and SAW 
are in different universality classes. Guttman and Bursill 
conjectured v = 4/5 in two dimensions (see however [^). 
This was proved by Majumdar |^ by relating the LERW 
problem to that of random spanning trees and showing 
that the fractal dimension of LERWs is the same as that 
of chemical paths in random spanning trees. In two 



spanning trees to the q ^ limit of the (/-state Potts 
model from conformal field theory 

Bounds on the expected number of erased steps have 
been obtained by Lawler |^ and it was found that in < 4 
dimensions the fraction of steps remaining unerased — > 
as the number of steps N —i- oo. Duplantier obtained 
exact exponents for the behaviour of the probability that 
k LERWs of length n starting at neighboring points do 
not intersect and also the winding number distribution of 
a LERW [§. Lawler has also shown that the LERW 
is equivalent to the Laplacian random walk model studied 
by Lyklema and Evertsz ||ll[]. 

Our interest in this paper is the distribution of sizes 
of erased loops for the LERW problem. The correspond- 
ing question for the totally random walks is the well- 
known Polya problem ||l^. For SAWs the problem is 
also the well-studied question of enumeration of poly- 
mer rings ijl^. A similar problem is encountered in the 
context of self-organised Eulerian walkers model and has 
been studied recently by Shcherbakov et al [|l^ . We show 
that the probability P{1) that an erased loop has perime- 
ter I, equals Pst{l)/l where Pst{l) is the probability that 
a loop of perimeter I is formed when a bond is added to 
a random spanning tree. For large I, P{1) and Pst{l) are 
expected to show power law behaviours, say P(l) ^ 
and Pst{l) ^ Then our result implies that 
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(1) 



We give scaling arguments to derive the exponent a in 
terms of the fractal dimension z of chemical paths on ran- 
dom spanning trees. For deterministic fractals, this ex- 
pression is modified, and involves also the ratio of Haus- 
dorff and spectral dimensions of the fractal. As a simple 
illustrative example we consider the Sierpinski Gasket, 
and calculate the exponents z and t directly from first 
principles. 

Consider a N step simple random walk lu — 
[{uj{0),uj{l), ...uj{N))], where uj{k) is the position of the 
random walker on the lattice after k steps. Let j be the 
smallest value such that uj{j) = oj{i) for < i < j < N. 
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We then obtain a new walk, Co = [(w(0), w(l), ■■w(j), + 
1), ...w{N))\, by deleting all steps between i and j. This 
process, corresponding to removing loops from lo in 
chronological order, is repeated till a j can no longer be 
found. The resulting walk is a LERW of length n < N . 

Consider a LERW, on a d-dimensional hypercubical L"^ 
torus, formed from a N step simple random walk by eras- 
ing loops. We define P(Z, N, L) to be the probability that 
the {N + l)th step results in erasing a loop of perimeter 
I. Let 

P{1) = lim lim Pil,N,L) (2) 

L — >oo N — >oo 

Consider the random walk starting at O. After N steps 
of the walk, we consider the directed tree formed using 
last exit bonds from all the sites visited by the walk, 
except the endpoint of the walk. This is called the last 
exit tree Tjv after TV steps. For N » L'^, all sites of 
the lattice are visited at least once, and is a spanning 
tree. It was proved by Broder that in the steady state 
all such spanning trees occur with equal probability. The 
LERW after N steps is just the directed path from O to 
the endpoint of the walk along Tjv. 

Now, consider a particular loop C of I directed bonds 
bi,b2...bi (see Fig.l). Let P{£, bj) be the probability that 
the {N + l)th step of the walk will result in formation 
of the loop C in the LERW problem with the {N + l)th 
step being along the bond bj . This occurs if and only if 

(i) the {N + l)th step forms the loop £ on T^r, with bj 
as the last step, 

(ii) there is a directed path in from O to the head 
site of bj, which does not include any bonds in C. 

Let Pst{C,bj) denote the probability that (i) occurs. 
This probability is easy to compute using the break- 
collapse method jl^ collapsing the loop £ to a single 
point. Hence it is easy to to see that Pst{C,bj) is the 
same for all j from 1 to L Thus 

Pst{C,b,)^Pst{C)/l, (3) 

where Pst{^) denotes the probability that loop £ is 
formed on T^r, whatever the position of the last step. To 
calculate P{£, bj), we have to multiply Pst{C-i bj) by the 
conditional probability P{0\C,bj) that (ii) occurs given 
that loop C is formed on the spanning tree with bj as the 
last bond. Thus we have 

P{C,b,) = P,^{C,b,).P{0\C,b,) (4) 

As for any spanning tree Tjv with end point of walk on 
£, the directed path from O must lead to one of the sites 
in the loop C, we must have 

I 

Y,P{0\C,b,)^\ . (5) 

Summing Eq.(4) over j from 1 to /,and using Eq.(5) we 
get 



f. 

Y,P{C,b,) = P,,{C)/t (6) 
i=i 

Finally, we sum over different shapes and positions of the 
loop C having the same perimeter I, to get 

m = Pst{i)/i (7) 

In deriving this result we have used the fact that 
iV — > oo limit is taken before the L ^ oo limit. It 
seems reasonable that the order of limits in the defini- 
tion (1) can be interchanged without affecting the value 
P{1). However, our proof uses the spanning tree prop- 
erty, and hence needs modification ii L ^ oo limit is 
taken before the N ^ oo limit. 

On a square lattice it is easy to calculate Pst{l) exactly 
for smah values of Z 0- We thus find P(2) = 0.25, 
P(4) « 0.03681 and P(6) 0.01034. We have done 
Monte-Carlo simulations and verified these figures to very 
high accuracy. In two dimensions, t — 8/5 [|l^, and this 
implies that a = 13/5 for d ^ 2. This is also in good 
agreement with our simulations. 

For d < 4, the number of steps of the random walk of 
N steps that are still not erased is a negligible fraction 
of Then the expected number of erased steps per 

step of the random walk is 

oc 

IP{1) = 1 for d < 4 (8) 

1=2 

For d > 4, there is a finite probability that a bond of 
the random walk which is generating the LERW will not 
be erased at any future time. Let this probability be 
called Poo- Then for d > 4, the average length of loop 
erased walk for a random walk of A^ steps increases as 
NPao- Since the average lengh erased is N^'^^IP{1), 
this implies that 

oc 

^?P(0 = 1-Poo. (9) 

1=2 

Thus, using our relation between the LERWs and span- 
ning trees, we are led to the interesting and paradoxical 
result that Poo can be thought of as the probability that 
adding a bond at random to a random spanning tree will 
not form a loop of finite perimeter, and this is nonzero if 
the dimension d of the space in which the spanning tree is 
embedded is greater than 4. This result clearly depends 
on the fact that the thermodynamic limit of large system 
size is taken before the limit / — > oo in the summation. 
For finite lattices, adding a bond to a spanning tree must 
lead to formation of a loop. 

We now present a scaling argument to determine a in 
terms of the fractal dimension of the LERW {z = \/v). 
Let n{l,N) = no. of loops of length / generated when 
random walk is of A^ steps. The typical excursion of a 
random walk of A^ steps R varies as N^^^. 
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For 1 < d < 4, the linear size of the largest loop gen- 
erated is expected to be of the order of R. Since I ~ i?^ 
the perimeter of the largest loop ^ N^^^. For large N, 
n{l,N) grows as NP{1). From finite size scaling theory, 
we expect that for large / and N, n{l,N) satisfies the 
scaling form 



n{l,N) 



(10) 



where f{x) is a scaling function. For the cumulative dis- 
tribution [no. of loops of size > I] 



N 



T9i 



I 



(11) 



The scaling function g(x) is assumed to be finite at small 
X and decay rapidly for large x i.e. for I larger than the 
cutoff length iV^/^. For I = kN^^-^, where A; is a finite 
constant of order 1 we must have c{l) of order 1. This 
gives 
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z 



(12) 



In one dimension, the scaling argument given above 
breaks down. On a linear chain, the erased loops can 
only be of size 1 — 2. On more complicated but linear 
graphs, such as ladders, we can have loops of arbitrarily 
large values of I, but P(/) decays exponentially with 
and the size of the largest loop generated does not scale 
as i?, as assumed in the scaling argument. For d > 4, the 
LERW is approximately a random walk and a = (i/2. 

Note that the scaling relation Eq.(12) does not involve 
the dimension of space d explicitly, but still is valid only 
for d less than the the upper critical dimension 4. It is 
interesting to ask how this relation needs to be modified 
to remain valid for noninteger values of d. We confine 
our arguments to recursively defined fractals, which are 
explicitly constructed spaces having a noninteger d. 

The scaling argument above is easily extended to work 
for deterministic fractals. In this case, for a random walk 
of iV steps, R iV'/^d (fQj. d < 2), where d and d are 
the spectral and hausdorff dimensions of the fractal. Re- 
peating the above argument, we then obtain 



a = 1 



2d 



for d<2. 



(13) 



Calculation of the chemical distance exponent z for 
spanning trees is quite straightforward for simple deter- 
ministic fractals of finite ramification index. Since this 
calculation has not appeared so far in the literature we 
describe it briefly below for the Sicrpinski Gasket (SG). 
The exact renormalization equations for spanning trees 
on the SG may be deduced from the general recursion 
equations for the q-state Potts model in the limit g — > 
p^ . However, for our purpose here, it is more convenient 
to use the recursion equations written down by Knezevic 



and Vannimenus (KV) iQ in the context of studying col- 
lapse transition of branched polymers on the SG. Only 3 
of the 6 graphs studied by KV have no vacant sites, and 
thus only these have non-zero weights for the problem of 
spanning trees (with no other interactions). These corre- 
spond to the cases where all the three vertices of the rth 
order triangle are connected to each other using bonds 
within the tree, two are connected to each other and not 
to the third, and all three are unconnected [Fig. 2] Let 
these weights be called A^'^\ B^'"-' and C*-'"-' respectively. 
By definition A^"^^ gives the number of spanning trees 
on the rth order gasket, B^^^ gives the number of two- 
rooted trees with two vertices connected and C^'"' gives 
the number of three-rooted spanning trees with all ver- 
tices unconnected. From KV, or directly, the recursion 
equations for A,B and C are easily written down 



(14) 



The initial values are given by A^^) = 3,B^^'> = l,C(i) = 
1. We define a new variable X^"^^ = A('')C(''V-B('')' . It 
is easy to see that X^^^ satisfies the following recursion 
equation: 



2XW + 7/3 



49/36 + XW736 + 7XM /18 



(15) 



This equation has the fixed point X* = 3. Let la and lb be 
the average lengths of the chemical paths connecting the 
lower two vertices of the rth order generating functions 
A and B respectively. To find the recursions for la and 
lb consider, for example, the graph of order (r -I- 1) and 
type B shown in Fig. 3 (a). The probability of this graph 
is ^WsMV-S^''^^^ = 1/(7 + X*) = 1/10. The length of 
the chemical path connecting the vertices is la+2lb. Thus 
the contribution of this graph to ll^~^^^ is {la + 24)/10. 
Summing up over all relevant graphs we get: 



^('■+1) = 5la/3 + 2lb/3 

t^'^ = + lb 



(16) 



We thus obtain la, Ih ^ X\ X = {20 + x/205)/15 being the 
largest eigenvalue in Eq.(16). Thus I ~ i?^ where R is 
the linear size of the gasket and z = In A/In 2. Putting in 
the value of A we get z = 1.1939.... 

We now find the exponent r. We show in Fig. 3. how a 
loop of the order of R^ may be obtained by adding a bond 
so that the one of the lower order B type graph becomes 
a A type graph. Supposing there are R^ positions where 
we could have added the bond in order to get the loop. 
Then the probability of this event ^ R^ / R'^. Hence we 
obtain: 



(17) 
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Thus we can find r if we can determine the exponent 
/3, which gives the fractal dimension of the boundary be- 
tween the two constituting sub-trees of the B type graph. 
We note that times B^^^ gives the number of ways of 
getting a A graph by addition of a bond to a i? graph 
with the added bond labeUed. But for every resulting A 
graph the labelled bond could be anywhere on the back- 
bone (path joining the three corner vertices) of length of 
order R^. This gives us 



(18) 



Now we note that the resistance between two points on 
a lattice with unit resistances on all bonds is given by 
the ratio of number of two-rooted spanning trees, with 
roots at the two given points, to the number of single- 
rooted spanning trees. It follows then that the ratio B/A 
gives the resistance between the corner points of a trian- 
gle, which scales as R" and so is related to the spectral 
dimension of the lattice. It can be shown easily that 
a = 2d/d- d. Thus from Eq.(18) we get 



(3 = z — a = z — 2d/d + d 
Using this in Eq.(17) we get: 

T = 2d/zd 



(19) 



(20) 



From Eq.(7) the LERW exponent a = t+1 and we verify 
the result Eq.(13) obtained by simple scaling arguments. 
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Figure Captions 

FIG. 1. A spanning tree Tjv with the endpoint of walk at 
X. The loop C (shown in bold) of 10 directed bonds is formed 
if the bond 6g, denoted by the white arrow, is added to Tat. 
In this case there is a directed path from the origin O to the 
head of 6g. 



FIG. 2. Diagrams representing the generating functions 
for spanning trees on the Sierpinski gasket. 



FIG. 3. Formation of a loop of perimeter of the order of 
on addition of a bond. 
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